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RESUME 


Ce  document  Etudie  1* application  a  la  compression  d* impulsion 
numErique  de  deux  types  de  codage  de  la  friquence:  les  codes  de 
Costas  et  les  codes  £  congruences  quadratiques.  Ces  codes  ont  la 
particularity  de  produire  des  diagrammes  d'ambigulte  en  forme  de 
punaise  avec  de  faibles  lobes  secondaires.  Ce  rapport  traite  des 
diffErents  algorithmes  de  construction,  de  1' Evaluation  des 
fonctions  d'auto-corrElation  et  de  correlation  croisEe,  et  de  leur 
comparaison  £  des  codes  connus.  Bien  qu'ayant  des  rEsultats 
intEressants  pour  le  nombre  de  codes  disponibles  et  le  niveau  des 
lobes  secondaires,  ils  semblent  Etre  moins  per formants  sur  ces 
points  que  les  codes  pseudo-alEatoires,  "chirp”  et  "step-chirp". 


ABSTRACT 


In  this  report  we  present  the  analysis  and  construction 
methods  for  two  types  of  frequency  coded  pulses:  (a)  Costas  codes 
and  (b)  quadratic  congruential  codes.  These  codes  have  the 
characteristic  of  producing  a  thumbtack  ambiguity  diagram  with  low 
sidelobes.  Computer  simulations  were  performed  to  calculate  their 
autocorrelation  functions.  Comparative  evaluation  of  the  perfor¬ 
mance  of  these  codes  and  well  known  codes  was  carried  out.  Results 
showed  that  Costas  codes  and  quadratic  congruential  codes  have  low 
sidelobes  but  not  as  low  as  the  chirp,  step-chirp  or  pseudo-random 
codes . 
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EXECUTIVE  SUMMARY 


La  compression  d* impulsion  numdrique  peut  etre  utilisde  dans 
un  radar  &  fonctions  multiples  comme  contre-contre-mesure  dlectro- 
nique  pour  contrer  le  brouillage.  La  selection  aldatoire  d'un 
nouveau  code  d  chaque  transmission  permet  de  detecter  les  retours 
possddant  le  mime  code  que  celui  de  1* impulsion  transmise  tout  en 
rejetant  les  retours  ayant  tout  autre  code.  II  est  important  de 
choisir  les  codes  pour  que  leur  fonction  d'auto-corrdlation  ait  un 
lobe  principal  etroit  et  des  lobes  secondaires  trds  faibles,  et  que 
leur  correlation  croisde  soit  beaucoup  plus  faible  que  leur  auto¬ 
correlation  respective. 

Ce  document  etudie  1' application  a  la  compression  d* impulsion 
numdrique  de  deux  types  de  codage  de  la  frequence;  les  codes  de 
Costas  et  les  codes  d  congruences  quadratiques.  Ces  codes  ont  la 
particularite  de  produire  des  diagrammes  d'ambiguitd  en  forme  de 
punaise  tout  en  ayant  de  faibles  lobes  secondaires.  Ce  rapport 
traite  des  diffdrentes  rndthodes  de  construction,  de  devaluation 
sur  ordinateur  des  fonctions  d'auto-corrdlation  et  de  correlation 
croisde,  et  de  leur  comparaison  d  des  codes  connus,  tels  les  codes 
biphases  pseudo-aldatoires  et  les  codes  modules  en  frequence  chirp 
et  step-chirp. 

Les  fonctions  d'auto-corrdlation  et  de  correlation  croisde  ont 
dtd  examinees  pour  des  codes  de  longueur  voisine  de  128,  limite 
fixde  par  l'appareillage  ddjd  existant  et  utilisde  lors  de  prece- 
dentes  recherches  sur  d'autres  codes.  Suite  d  1* analyse  des 
rdsultats,  les  codes  de  Costas  ont  dtd  trouves  les  meilleurs  en  ce 
qui  concerne  la  fonction  d'auto-corrdlation,  surtout  avec  la 
mdthode  de  construction  de  Welch,  alors  que  les  codes  quadratiques 
offrent  de  meilleurs  rdsultats  pour  la  correlation  croisde. 
Cependant,  les  fonctions  d'auto-corrdlation  de  ces  codes  possddent 
des  lobes  secondaires  plus  dlevds  que  celles  des  codes  chirp,  step- 
chirp  et  pseudo-aldatoires  equivalents. 

Les  codes  de  Costas  et  les  codes  quadratiques  peuvent  etre 
utilises  pour  la  compression  d' impulsion  numdrique  mais  ne  sont  pas 
aussi  performants  que  les  codes  pseudo-aldatoires,  chirp  et  step- 
chirp  en  ce  qui  a  rapport  au  niveau  des  lobes  secondaires  des 
fonctions  d'auto-correiation  et  au  nombre  de  codes  disponibles. 
Les  codes  quadratiques  sont  plus  approprids  pour  une  utilisation 
ndcessitant  un  compromis  entre  une  bonne  fonction  d ' auto-corre¬ 
lation  et  une  bonne  correlation  croisde. 
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1.  INTRODUCTION 


Les  signaux  codes  sont  utilises  pour  les  radars  a  compression 
d' impulsion.  Ils  permettent  d'am61iorer  la  resolution  et  la  detec¬ 
tion  des  cibles  sans  a  avoir  augmenter  la  puissance  crete  emise. 
Les  impulsions  transmises  sont  modulees  en  frequence  ou  en  phase 
selon  un  code  specif ique.  A  leur  retour,  les  impulsions  sont  trai- 
t6es  avec  un  filtre  adapt e  au  code  utilise  lors  de  la  transmission. 
La  sortie  du  filtre  represente  la  fonction  de  correlation  entre 
1' impulsion  transmise  et  1' impulsion  regue.  Pour  pouvoir  detecter 
des  signaux  faibles  proches  de  signaux  plus  forts,  la  fonction 
d# auto-correlation  du  code  devrait  idealement  comporter  un  pic 
principal  dleve  et  etroit,  et  des  pics  secondaires  faibles  et  tous 
de  mAme  amplitude.  En  pratique,  on  obtient  souvent  des  pics 
secondaires  elev£s  que  l'on  doit  chercher  a  minimiser. 

Dans  certaines  applications,  le  radar  utilise  l'agilite  de 
codage  qui  consiste  a  utiliser  un  code  different  a  chaque  transmis¬ 
sion.  Cette  technique  de  contre-contre-mesure  electronique  permet 
de  detecter  les  cibles  possedant  le  meme  code  que  celui  de 
1' impulsion  transmise,  et  de  rejeter  le  brouillage  et  les  echos 
lointains  dont  les  cibles  ont  un  code  different.  II  est  alors 
necessaire  que  la  correlation  crois&e  entre  les  codes  successifs  ne 
comporte  que  des  pics  faibles. 

La  correlation  croisde  pourra  admettre  des  pics  dominants  a 
condition  qu'ils  soient  decentres  pour  ne  pas  ressembler  a  une 
auto-correlation.  D/une  impulsion  a  1# autre,  les  pics  principaux 
des  cibles  renvoyant  le  bon  code  sont  tou jours  au  meme  endroit, 
tandis  que  les  pics  dominants  des  cibles  possedant  un  code 
different  sont  detectees  a  des  positions  diverses.  La  detection 
cumulative  ou  1 ' integration  coh6rente  permet  de  concentrer  la 
detection  des  retours  avec  le  bon  code,  tout  en  eparpillant  les 
detections  des  retours  dont  le  code  ne  correspond  pas  a  celui 
transmis  [1,2,3]. 

Il  existe  diff6rentes  methodes  de  codage  de  signaux.  Une 
modulation  discrete  de  la  frequence  peut  etre  obtenue  en  6chantil- 
lonnant  une  modulation  continue  de  la  frequence.  Egalement, 
l'6chantillonnage  de  la  phase  d'un  signal  module  en  frequence  peut 
etre  utilise  pour  un  codage  de  la  phase.  Nous  etudierons  quelques 
types  de  codage  en  frequence  suggeres  pour  la  compression  d' impul¬ 
sion  tels  que  les  codes  de  Costas  [4, 5, 6, 7]  et  les  codes  a 
congruences  quadrat iques  [8,9,10]  et  les  comparerons  a  d'autres 
codes  deja  utilises  en  radar  comme  les  codes  chirp,  step-chirp  et 
pseudo-aieatoires  [11,12,1].  Tous  les  codes  etudies  ont  une 
longueur  voisine  de  128  car  1 'appareillage  disponible,  construit 
lors  d' etudes  pr6cedentes  sur  d'autres  codes,  peut  produire  des 
codes  ayant  au  maximum  128  points  [12,13], 
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Pour  etudier  les  proprieties  des  codes  et  determiner  quel  type 
donne  les  meilleures  fonctions  d' auto-correlation  et  de  correlation 
croisde,  nous  avons  construit  un  logiciel  qui  simule  la  correlation 
entre  les  signaux  transmis  et  regu.  En  premier,  un  algorithme 
different  pour  chaque  code  est  utilise  pour  calculer  la  frequence 
et  la  phase  instantanees  du  signal  desire.  Par  la  suite,  le 
programme  simule  l'echantillonnage  des  composantes  en  phase  et  en 
quadrature  (I  &  Q)  d'un  signal  en  produisant  une  sequence  de 
variables  complexes.  Finalement,  des  transformees  de  Fourier 
rapide  sont  employees  pour  effectuer  la  correlation  entre  les 
sequences  transmise  et  regue.  Le  logiciel  calcule  egalement  la 
valeur  du  pic  central,  du  pic  secondaire  maximum  et  la  moyenne  des 
pics  secondaires  de  la  correlation.  Pour  chaque  fonction  d' auto¬ 
correlation  des  codes  etudies,  la  largeur  et  la  hauteur  du  pic 
principal  sont  tou jours  respect ivement  egales  a  un  et  a  la  longueur 
du  code. 
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2.  LES  CODES  DE  COSTAS 


•  Les  codes  etudies  dans  ce  rapport  s/appliquent  h  une  modula¬ 

tion  discrete  de  la  frequence.  Les  impulsions  sont  divisees  en 
segments  et  la  frequence  du  signal  varie  avec  chacun  des  segments. 

A  Une  matrice  est  utilisee  pour  representer  le  codage  en  frequence 

d'un  signal.  Si  a(n)  represente  la  position  de  1* element  dans  la 
matrice  du  code,  la  frequence  de  chaque  segment  est  donnee  par: 

f (n)  =  fQ  +  a(n)  •  B/N  (  1  ) 

ou  f0  est  la  frequence  initials  du  signal,  B  la  largueur  de  bande 
du  signal,  et  N  le  nombre  de  segments. 

Un  code  de  Costas  est  une  matrice  representant  la  frequence  en 
fonction  du  temps  pour  le  signal  code.  La  matrice  comme  celle  de 
la  figure  1  ne  possede  qu'une  seule  frequence  pour  chaque  inter- 
valle  de  temps  et  un  seul  intervalle  de  temps  pour  chaque  fre¬ 
quence.  II  y  a  N  elements  presents  dans  une  matrice  d'un  code  de 
Costas  et  la  position  de  ces  elements  est  une  permutation  de 
l'ensemble  {l,2,...,N}. 


j 


Figure  1  Code  de  Costas  (m6thode  de  Welch?  q=ll,a=6, rj=2 )  . 

Un  code  repond  au  critAre  de  Costas  si,  en  superposant  la 
matrice  et  une  version  decaiee  d'elle-meme,  au  plus  deux  elements 
coincident,  et  ce  pour  n'importe  quel  decalage  non  nul  en  temps 
et/ou  en  frequence.  Le  critere  de  Costas  permet  aux  codes  de  Welch 
et  de  Golomb,  ainsi  qu'aux  autres  codes  de  Costas,  d' avoir  un 
diagramme  d'ambiguite  et  une  fonction  d'auto-correiation  quasi- 
ideaux,  c'est-A-dire  une  pointe  centrale  et  des  lobes  secondaires 
d'amplitude  faible  et  constante  sur  toute  la  region. 

Les  constructions  de  codes  de  Costas  etudiees  sont  basees  sur 
la  theorie  des  corps  de  Galois  [14].  Un  corps  de  Galois  est  un 
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ensemble  dont  la  dimension  (q)  est  un  nombre  premier  tel  qu'il 
existe  l'ensemble  GF(q)  *  {0,1,2, ... ,q-l)  sur  lequel  les  operations 
arithmdtiques  se  font  modulo  q.  Un  dldment  primitif  d'un  corps  de 
Galois  permet  d/ exprimer  tous  les  dldments  non  nuls  de  1' ensemble 
par  une  puissance  de  ce  nombre (modulo  q) .  Chaque  ensemble  fini 
possede  plusieurs  elements  primitif s.  Soit  un  corps  de  Galois 
GF(q)  oil  q,  la  dimension  de  1  •  ensemble ,  est  un  nombre  premier.  On 
peut  decomposer  le  nombre  q-1  par  une  factorisation  de  nombre 
premier: 


(q-1)  =  Pi811  •  p/*  .....  pr~  (  2  ) 

Le  nombre  d'dldments  primitifs  de  l'ensemble  GF(q)  est  alors  donnd 
par  : 

4>  (q-1)  =  Pi"1'1  •  Pi02'1  •  .  •  .  •  Pr”'1  •  (Pi  -  1)  •  •  •  •  •  (Pr  -  1)  (  3  ) 

Par  exemple,  pour  un  ensemble  GF(ll),  le  corps  de  Galois 
servant  &  generer  le  code  est  1' ensemble  suivant: 

GF(ll)  «  {0,1,2, 3, . . . ,10} 

Dans  ce  cas-ci,  q=*ll;  on  peut  done  exprimer  par  factorisation 

(q-1)  =  10  =  21  •  51 

Le  nombre  d'616ments  primitifs  de  ce  corps  de  Galois  se  calcule 
comme  suit  d'apr6s  1' Equation  (3): 

^ (q-1)  =  *(10)  -  2°  •  5°  • (1) • (4)  =  4 

L'ensemble  GF(ll)  possede  done  quatre  dldments  primitifs.  Le 
nombre  six  est  un  des  elements  primitifs  de  GF(ll)  car  les 
puissances  de  ce  nombre,  modulo  11,  redonnent  tous  les  elements  non 
nuls  de  l'ensemble  GF(ll)  =  {0,1,2, ... ,10} . 


61  =  6 
6Z  —  3 
63  *  7 
6*  *  9 
63  =  10 
6*  -  5 
67  -  8 

68  -  4 

69  *  2 
610  =  1 

Prenons  un  dldment  primitif  a  d'un  corps  de  Galois  GF(q) .  On 
remarque  que 

a<rl  |  mod  q  =  1 
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et  que  pour  n'importe  quel  entier  i, 

a  |  mod  q  =  (®  I  ®od  q  *  a)  |  mod  q 

ce  qui  permet  de  calculer  les  puissances  d'un  Element  primitif  par 
recurs ivitd  [15].  Cette  mgthode  pr6sente  1 'avantage  de  tou jours 
traiter  des  nombres  plus  petits  que  <*3. 

La  liste  des  616ments  primitifs  des  corps  de  Galois  utilises 
dans  ce  texte  est  donnde  dans  le  tableau  1. 


131  I  2,  6,  8,  10,  14,  17,  22,  23,  26,  29, 


30, 

31,  37 

,  40, 

50, 

54, 

56,  57 

,  66, 

67, 

72,  76 

,  82, 

83, 

85, 

87,  88 

,  90, 

93, 

95,  96 

,  97, 

98, 

103, 

,  104, 

106, 

110, 

r  111, 

115, 

116, 

118, 

,  119, 

120, 

122, 

r  124, 

126, 

127, 

128 

Tableau  1  Elements  primitifs  de  certains  corps  de  Galois. 

Diffdrents  types  de  constructions  permettent  d'obtenir  des 
codes  r£pondant  au  critere  de  Costas.  Nous  6tudierons  les  m£thodes 
de  Welch  et  de  Golomb.  Ces  codes  de  Costas  sont  g£n4r&s  h  partir 
d'un  corps  de  Galois  GF(q)  ou  (q-1)  et  (q-2)  representent  respec- 
tivement  le  nombre  d'616ments  pour  les  codes  de  Welch  et  de  Golomb. 


2.1  MfeFHODE  DE  WELCH 

D'apr&s  Welch,  on  obtient  un  code  de  Costas  si  a  est  un 
616ment  primitif  de  GF(q) ,  t\  est  un  616ment  non  nul  de  1 'ensemble 
et  q  est  un  nombre  impair  et  premier.  La  matrice  du  code  se 
construit  comme  suit: 


(i,j)  -  (i»  n  *  <*1 1  mod  q)  (  4  ) 

i  :  position  horizontale 
j  :  position  verticale 
1  s  i,j  s  q-1 

Cette  m6thode  produit  (q-1)  •  4>  (q-1)  codes  dif f erents  de 
dimension  (q-1) • (q-1) .  C'est-fc-dire  que  l'on  peut  g£n£rer  (q-1) 
codes  avec  le  m&me  616ment  primitif  en  variant  la  valeur  du 
coefficient  >7.  Ces  dif f 6rents  codes  g£n£ry  par  le  meme  Element 
primitif  ne  sont,  en  fait,  que  des  rotations  cycliques  l'un  de 
1 ' autre . 

Pour  illustrer  un  code  de  Welch,  prenons  le  corps  de  Galois 
GF(ll) .  La  longueur  du  code  g£n6r6  sera  de  10.  En  choisissant  6 
comme  616ment  primitif  et  un  coefficient  6gal  h  2,  nous  obtenons  la 
matrice  de  la  figure  1,  oh  les  valeurs  de  i  et  j  sont: 


i 

j-  2 

•61  bo* 

1 

2»6l 

*  1 

2 

2  •  62 

«  6 

3 

2»63 

■  3 

4 

2*6* 

»  7 

5 

2»  63 

»  9 

6 

2*6“ 

=  10 

7 

2  •  67 

«  5 

8 

2*6® 

»  8 

9 

2  •  6® 

a  4 

10 

2*610 

-  2 

Prenons  un  code  de  Welch  de  longueur  126  ayant  un  produit 
temps-largeur  de  bande  (BT)  de  126,  126  616ments,  le  nombre  6  comme 
£16ment  primitif  et  un  coefficient  de  22.  Pour  ce  code,  nous 
obtenons  des  pics  secondaires  de  1' autocorrelation  &  -24.0  dB  par 
rapport  au  pic  central.  Changeons  cette  fois  pour  un  element  de 
116  et  un  coefficient  de  111,  nous  obtenons  alors  un  pic  secondaire 
maximal  de  -20.3  dB  sous  la  valeur  du  pic  principal. 

II  y  a  au  total  4536  (36»126)  possibility  de  codes  de  Welch 
poss£dant  126  elements.  Pour  chacune  de  ces  possibility,  le  pic 
secondaire  maximum  se  situe  entre  -20.3  dB  et  -24.0  dB  du  pic 
central.  Le  graphique  de  la  figure  2  repr6sente  la  fonction 
d' auto-corr61ation  d'un  code  de  Welch  dont  les  paramhtres  longueur. 
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produit  BT  et  nombre  d' elements  sont  egaux  a  126,  1 'element 

primitif  du  code  est  6  et  le  coefficient  vaut  22. 


Figure  2  Auto-corr61ation  d'un  code  de  Welch. 

On  peut  avoir  d ' autres  codes  de  Welch  de  meme  longueur  mais 
avec  un  nombre  diffdrent  d' elements.  II  s'agit  d'utiliser  moins 
d' elements  et  de  prendre  plus  d'echantillons  par  element  pour 
conserver  le  meme  nombre  total  d'echantillons.  Cependant,  la 
longueur  du  code  doit  etre  un  multiple  du  nombre  d' elements,  et  le 
produit  BT  doit  rester  egal  a  la  longueur  du  code.  Nous  allons 
etudier  deux  autres  cas  de  codes  de  Welch. 

Premi Arement ,  un  code  de  longueur  et  de  produit  BT  de  126  avec 
18  AlAments  fourni  A  partir  d'un  corps  de  Galois  GF(19).  Lorsque 
l'616ment  primitif  est  2  et  le  coefficient  6gal  A  9,  le  pic 
secondaire  principal  est  a  -18.1  dB  par  rapport  au  pic  central. 
Pour  un  dldment  primitif  de  2  et  coefficient  de  11,  le  pic  est  a 
-13.8  dB.  Pour  les  108  possibilitds  de  codes  de  Welch  de  longueur 
et  produit  BT  126  possAdant  18  Elements,  la  valeur  des  pics  secon- 
daires  principaux  se  situe  entre  -13.8  dB  et  -18.1  dB  sous  le  pic 
central  de  la  fonction  d'auto-corrAlation. 

Dans  le  cas  de  codes  de  longueur  et  produit  BT  egaux  a  128  et 
un  nombre  d'AlAments  6gal  A  16  (fourni  &  partir  d'un  corps  de 
Galois  GF ( 17 ) ) ,  la  valeur  des  128  diffArents  codes  pouvant  etre 
crAAs  varie  entre  -15.0  dB  et  -20.0  dB.  Les  limites  sont  respec- 
tivement  obtenues  avec  des  AlAments  primitif s  de  6  et  5,  et  avec 
des  coefficients  de  13  et  7. 

II  est  aussi  possible  d'obtenir  des  codes  de  Costas  de 
longueur  infArieure  a  (q-1)  a  partir  d'un  corps  de  Galois  GF(q) 
[7].  II  s'agit  de  retrancher  des  elements  de  la  matrice  dont  la 
valeur  aquivaut  a  un  coin  de  la  matrice.  La  matrice  de  dimension 
10  correspondent  au  code  de  Welch;  q=ll,  a=6  et  r)= 2,  et  representee 
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k  la  figure  1,  peut  servir  d'exemple.  Le  premier  element  de  la 
matrice  vaut  1  et  le  dernier  vaut  2.  On  peut  alors  enlever  les 
deux  rangees  et  les  deux  colonnes  contenant  ces  deux  elements  pour 
obtenir  un  matrice  de  dimension  8. 

II  est  tou jours  possible  d' obtenir  un  code  de  Welch  de 
longueur  (q-3)  k  partir  du  corps  de  Galois  GF(q) ,  si  celui-ci  a  2 
comme  element  primitif  et  1  comme  coefficient.  Dans  certains  cas, 
elle  peut  etre  la  seule  possibility.  C'est  ainsi  que  nous  avons 
obtenu  un  code  de  Welch  de  longueur  128  k  partir  du  corps  de  Galois 
GF (131) . 

Nous  avons  etudie  les  propriet£s  de  ce  nouveau  code  en 
calculant  sa  fonction  d' auto-correlation.  Lorsque  le  produit  BT 
est  egal  k  la  longueur  du  code  et  le  nombre  d'elements,  soit  128, 
le  pic  secondaire  maximum  est  a  -21.5  dB  par  rapport  au  pic 
central.  Ces  valeurs  se  situent  dans  la  moyenne  des  resultats  des 
codes  de  Welch  precedemment  etudies. 

Le  tableau  2  contient  certaines  valeurs  de  la  fonction  d' auto¬ 
correlation  de  differents  codes  de  Welch  etudies  . 


longueur 

nombre 
d' elements 

meilleur 

a  r/  pic 

[  dB  ] 

a 

pire 

r)  pic 

[  dB  ] 

126 

18 

2 

9 

-18.1 

2 

11 

-13.8 

126 

126 

6 

22 

-24.0 

116 

111 

-20.3 

128 

16 

5 

7 

-20.0 

6 

13 

-15.0 

128 

_ 

128 

2 

1 

-21.5 

— 

—  — 

— 

Tableau  2  Auto-correlation  des  codes  de  Welch. 

Le  maximum  de  la  correlation  croisee  entre  deux  codes  de  Welch 
a  une  valeur  moyenne  de  -14  dB  par  rapport  au  pic  principal  de  la 
fonction  d' autocorrelation.  La  correlation  croisee  entre  deux 
codes  de  Welch  generes  par  le  meme  element  primitif  presente 
cependant  un  ou  deux  pics  importants  mais  decentres.  Ceci  est  du 
k  la  relation  entre  les  codes  ayant  le  meme  element  primitif  mais 
non  le  m£me  coefficient.  Le  graphique  de  la  figure  3  presente  la 
fonction  de  correlation  croisee  entre  deux  codes  de  Costas  generes 
par  le  m£me  element  primitif,  c'est-e-dire  deux  codes  de  Welch  de 
longueur  et  produit  BT  de  126  et  possedant  18  elements.  Les  deux 
sont  generes  par  un  element  primitif  de  2;  un  poss£de  un  coef¬ 
ficient  de  9  et  1' autre  un  coefficient  de  11. 


8 


welch  -  n:  126  ;bt:126  ;step:18  ;olpha:2  ;coef:  9 
welch  -  n:126  :bt:126  ;step:18  ;dpha:2  ;coef:  11 


Figure  3  Correlation  croisee  entre  deux  codes  de  Welch  generes  par 
le  meme  element  primitif. 

2.2  METHODE  DE  GOLOMB 

Avec  la  methode  de  Golomb,  on  peut  obtenir  O (q-1) ] 2  codes 
differents  de  dimension  (q-2) • (q-2).  Un  element  est  place  dans  la 
matrice  du  code  en  position  (i,j)  si  et  seulement  si: 

(a1  +  0J)  Imod,  -  1  (  5  ) 

i  :  position  horizontale 
j  :  position  verticale 
1  <  i,j  s  q-2 

Les  operations  se  faisant  modulo  q.  Les  parametres  q  et  jS 
sont  des  elements  primitif s  de  GF(q) ,  q  etant  un  nombre  impair  et 
premier. 


j 


123456789  i 


Figure  4  Code  de  Costas  (methode  de  Golomb;  q=ll,a=6,£=7) . 
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Avec  deux  elements  primitifs  du  corps  de  Galois  GF(ll)  tels 
que  a=6  et  £=7,  on  obtient  un  code  de  Golomb  de  longueur  9 
repr6sente  par  une  matrice  comme  celle  de  la  figure  4,  ou  les 
valeurs  de  i  et  j  sont: 


i 

j 

(6 

1  +  73) 

nod 

li 

= 

1 

1 

7 

61 

+ 

77 

SS 

6 

+ 

6 

= 

12 

= 

1 

2 

8 

62 

+ 

7® 

S 

3 

+ 

9 

= 

12 

= 

1 

3 

2 

63 

+ 

72 

= 

7 

+ 

5 

= 

12 

= 

1 

4 

4 

64 

+ 

74 

= 

9 

+ 

3 

= 

12 

= 

1 

5 

3 

63 

+ 

73 

10 

+ 

2 

= 

12 

= 

1 

6 

1 

6® 

+ 

71 

5 

+ 

7 

= 

12 

= 

1 

7 

6 

67 

+ 

7® 

= 

8 

+ 

4 

= 

12 

= 

1 

8 

9 

6® 

+ 

7® 

= 

4 

+ 

8 

= 

12 

= 

1 

9 

5 

6® 

+ 

7s 

2 

+ 

10 

= 

12 

= 

1 

La  methode  de  construction  de  Golomb  genhre  des  codes  de 
Costas  diffdrents  de  ceux  produits  par  la  methode  de  Welch. 
Prenons,  par  exemple,  un  code  de  Golomb  fourni  a  partir  d ' un  corps 
de  Galois  GF(127)  oh  les  paramfetres  longueur,  produit  BT  et  nombre 
d' elements  sont  egaux  4  125.  Le  maximum  des  pics  secondaires  des 
1296  possibilites  de  codes  varie  alors  entre  -12.9  dB  pour  des 
616ments  primitifs  de  57  et  106,  et  de  -18.6  dB  pour  des  elements 
primitifs  de  7  et  55. 

Le  nombre  d'41£ments  des  codes  de  Golomb  peut  etre  different 
de  la  longueur  du  code.  Pour  les  16  codes  de  Golomb  de  longueur  et 
de  produit  BT  126  ayant  9  elements  (fourni  4  partir  d'un  corps  de 
Galois  GF ( 11) ) ,  le  maximum  des  pics  secondaires  est  de  -8.9  dB  sous 
le  pic  central,  et  ce  peu  importe  la  valeur  des  elements  primitifs 
du  code.  Par  contre,  la  valeur  moyenne  des  pics  secondaires  varie 
entre  -27.5  dB  et  -27.8  dB  respectivement  pour  des  elements 
primitifs  de  6  et  6,  et  de  2  et  8.  Les  codes  precedemment  cites 
sont  cependant  les  seules  possibilites  de  codes  de  Golomb  avec  un 
nombre  d'414ments  non  6gal  4  la  longueur  du  code. 

Le  mdme  principe  que  celui  utilise  pour  les  codes  de  Welch 
peut  4tre  appliqu6  pour  obtenir  des  codes  de  Golomb  de  longueur 
inf6rieure  4  (q-2) .  Si  a  +  ^  =  1  (a  et  £  6tant  deux  elements 
primitifs  de  GF (q) ) ,  le  premier  61ement  de  la  matrice  vaut  1;  la 
premiere  rangee  et  la  premiere  colonne  peuvent  alors  etre  enlevees 
pour  g6n6rer  un  code  de  Golomb  de  longueur  (q-3)  au  lieu  de  (q-2) . 

Cette  m6thode  permet  de  construire  12  codes  differents  de 
Golomb  de  longueur  128  4  partir  de  GF(131) .  Les  r£sultats  obtenus 
pour  la  fonction  d'auto-corr£lation  de  ces  codes  sont  semblables  4 
ceux  des  autres  codes  de  Golomb  d£j4  etudies.  Le  pic  secondaire 
maximum  se  situe  ici  entre  -13.8  dB  et  -17.7  dB  par  rapport  au  pic 
central  pour  un  produit  BT  6gal  4  la  longueur  du  code (128)  et  au 
nombre  d'616ments. 
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Les  rdsultats  obtenus  pour  1' auto-correlation  des  codes  de 
Golomb  avec  notre  algor ithme,  sont  presentes  au  tableau  3. 


longueur 

nombre 
d  Elements 

a 

minimum 

0  pic 

C  dB  ] 

a 

maximum 

0  pic 

[  dB  ] 

125 

125 

7 

55 

-18.6 

57 

106 

mmm 

126 

9 

tous 

tous 

-  8.9 

— 

— 

128 

128 

17 

115 

-17.7 

26 

106 

Tableau  3  Auto-correlation  des  codes  de  Golomb. 


V 
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3.  LES  CODES  A  CONGRUENCES  QUADRAT I QUES 


Les  codes  de  Welch  et  de  Golomb  sont  des  codes  de  Costas  car 
ils  repondent  au  critere  de  Costas.  Un  type  different  de  construc¬ 
tion  de  codes  pour  la  modulation  discrete  de  la  frequence,  se  base 
sur  la  theorie  des  congruences  quadratiques ;  ce  sont  ce  que  nous 
appellerons  les  codes  quadratiques.  Contrairement  aux  codes  de 
Costas,  la  sequence  d'un  code  quadrat ique  de  longueur  N  n'est  pas 
une  permutation  de  1' ensemble  {1,2, ...,N}.  Ces  codes  sont  de  type 
non-complet,  certaines  frequences  sont  presentes  plus  d'une  fois  et 
d'autres  ne  sont  jamais  utilisAes. 

Une  matrice  represente  les  valeurs  discretes  de  la  frequence 
du  code.  La  position  de  chaque  element  (i,j)  de  la  matrice 
representant  un  code  quadrat ique  de  longueur  N  est  donnee  par: 

j  *  I"  &±h.Li.-LLY  1  (  6  ) 

L  2  J  5,0,)  N 

1  s  i  s  N 
1  s  a  <  N-l 
a  :  coefficient 

Avec  cette  methode,  on  peut  obtenir  (N-l)  codes  differents  de 
longueur  N.  Les  matrices  des  codes  quadratiques  ont  la  propriety 
d'etre  symetriques  par  rapport  A  leur  centre,  c'est-A-dire  que 

(i,j)  =  (N+l-i,j)  pour  tout  N  impair 

Pour  bien  voir  la  symetrie  des  codes  quadratiques,  regardons 
a  la  figure  5  la  matrice  d'un  code  de  longueur  9  ayant  un  coef¬ 
ficient  egal  A  5. 


j 


Figure  5  Code  quadratique;  N=9 ,  a=5 . 
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oil  les  valeurs  de  i  et  j  sont: 


V 


i  j=  5»  i*  (i-l)/2  |  mod  9 


1 

2 

3 

4 

5 

6 

7 

8 
9 


5- 1*0/2 
5* 2* 1/2 
5* 3 *2/2 
5*4*  3/2 
5* 5* 4/2 
5*6*  5/2 
5* 7* 6/2 
5*8*7/2 
5 *9* 8/2 


0  =  0 
5  =  5 
15  =  6 
30  =  3 
50  =  5 
75  =  3 
105  =  6 
140  =  5 
180  =  0 


Puisque  la  longueur  du  code  doit  etre  un  nombre  impair  et 
premier,  nous  avons  etudie  les  codes  quadratiques  de  longueur  127. 
Les  meilleurs  resultats,  c'est-&-dire  les  pics  secondaires  les  plus 
faibles,  ont  dte  obtenus  pour  un  coefficient  de  valeur  egale  a  15. 
Le  pic  secondaire  principal  est  alors  de  -18.0  dB  par  rapport  au 
pic  central,  et  ceci  pour  un  produit  BT  egal  &  la  longueur  du  code 
ainsi  qu'au  nombre  d' elements,  soit  127.  Par  contre,  pour  un 
coefficient  de  107,  le  pic  secondaire  est  &  -15.8  dB  par  rapport  au 
r>ic  central.  Pour  toutes  les  autres  valeurs  du  coefficient,  entre 
1  et  127,  le  niveau  maximum  des  pics  secondaires  se  situe  entre  ces 
deux  bornes. 


Dans  1' article  «Time-frequency  hop  signals  part  II:  coding 
based  upon  quadratic  congruences»  [10]  ,  1' auteur  utilise  aussi  des 
nombres  N  impairs  mais  non  premiers  pour  gdnerer  des  codes  quadra¬ 
tiques  de  longueur  N.  Etant  donne  la  propriety  qu'ont  les  codes 
quadratiques  d'etre  symetriques  par  rapport  a  leur  centre,  nous 
g6n6rons  un  code  de  longueur  128  &  partir  d'un  code  de  longueur  129 
(impair,  non  premier)  duquel  nous  retranchons  1 'element  central 
pour  obtenir  128  elements.  Nous  avons  ensuite  etudie  les  proprie¬ 
ty  de  ces  codes  quadratiques  pairs.  Lorsque  le  produit  BT  et  le 
nombre  d'€l£ments  sont  tous  deux  dgaux  &  la  longueur  du  code,  ici 
128,  la  valeur  pointe  des  pics  secondaires  varie  entre  -0.4  dB  et 
-14.2  dB  sous  le  pic  central. 


Nous  avons  aussi  dtudie  quelques  autres  cas  de  codes  quadrati¬ 
ques  pairs  dont  le  nombre  d' elements  est  un  facteur  de  la  longueur 
du  code.  Nous  obtenons  pour  un  code  quadrat ique  de  longueur  et  de 
produit  BT  126  avec  18  elements,  des  pics  secondaires  principaux 
entre  -9.9  dB  et  -14.1  dB  par  rapport  au  pic  central.  Si  le  nombre 
d'616ments  est  dgal  h  14,  les  bornes  des  pics  secondaires  prin¬ 
cipaux  sont  -2.2  dB  et  -8.4  dB.  Prenons  des  codes  quadratiques  de 
longueur  et  de  produit  BT  de  128,  avec  success ivement  8  et  16 
616ments,  on  obtient  respect ivement  des  limites  de  -1.5  dB  et  -5.0 
dB  et  de  -9.8  dB  4  -13.8  dB  pour  les  maximums  des  pics  secondaires 
principaux. 
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Les  codes  guadratiques  etudies  possedent  entre  eux  une 
fonction  de  correlation  croisee  dont  la  valeur  moyenne  du  maximum 
est  h  un  niveau  de  -15.0  dB  par  rapport  au  pic  central  de  leur 
fonction  d' auto-correlation.  Les  meilleurs  resultats  sont  par 
contre  obtenus  avec  les  codes  guadratiques  pairs  de  longueur, 
produit  BT  et  nombre  d' elements  egaux  a  128.  La  fonction  de 
correlation  croisee  pour  ces  codes  est  a  -23.4  dB  par  rapport  au 
pic  principal  de  leur  fonction  d' auto-correlation.  Un  exemple  de 
fonction  de  correlation  croisee  entre  deux  codes  guadratiques  pairs 
de  longueur,  produit  BT  et  nombre  d' elements  128,  possedant  des 
coefficients  de  41  et  86  est  illustree  sur  le  graphique  de  la 
figure  6. 


Figure  6  Correlation  croisee  de  codes  quadratiques. 

Le  tableau  4  donne  quelques  resultats  pour  des  codes  possedant 
un  produit  BT  egal  a  leur  longueur.  Les  codes  dont  le  nombre 
d' elements  est  un  facteur  de  la  longueur  du  code  n' off rent  pas  de 
meilleurs  resultats  ni  au  point  de  vue  de  la  correlation  croisee, 
ni  pour  1 'auto-corrdlation. 


longueur 

nombre  d' 
Elements 

minimum 
t)  pic 

[  dB  ] 

maximum 
n  pic 

[  dB  ] 

126 

14 

1 

-  8.4 

5 

-  2.2 

126 

18 

5 

-14.1 

15 

-  9.9 

126 

126 

13 

-16.1 

88 

-12.5 

127 

127 

15 

-18.0 

107 

-15.6 

128 

8 

2 

-  5.0 

6 

-  1.5 

128 

16 

2 

-13.8 

10 

-  9.8 

128 

128 

41 

-14.2 

86 

-  0.4 

Tableau  4  Autocorrelation  des  codes  quadratiques. 
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4.  COMPARAI SON  DES  DIFFERENTS  TYPES  OE  CODES 


Nous  avons  examine  les  caracteristiques  des  trois  differents 
types  de  codes.  Les  codes  de  Costas  ont  une  fonction  d# auto¬ 
correlation  offrant  de  meilleurs  rdsultats  que  celle  des  codes 
quadratiques.  Parmi  les  codes  de  Costas,  les  resultats  obtenus 
montrent  que  les  codes  de  Golomb  offrent  de  moins  bonnes  possibi¬ 
lity  que  les  codes  de  Welch  pour  la  fonction  d' auto-correlation. 
Les  meilleurs  rdsultats  obtenus  par  les  codes  de  Golomb  sont  moins 
bons  que  n'importe  quel  resultat  obtenu  avec  les  codes  de  Welch 
lorsque  le  produit  BT  est  egal  a  la  longueur  du  code.  Par  exemple, 
pour  des  codes  possedant  une  longueur,  un  produit  BT  et  un  nombre 
d '  elements  egaux  a  128,  les  resultats  pour  un  code  quadrat ique  sont 
de  -13.8  dB,  de  -17.7  dB  pour  un  code  de  Golomb  et  de  -21.5  dB  pour 
un  code  de  Welch. 

En  ce  qui  concerne  la  correlation  croisee,  les  resultats  sont 
cependant  semblables  pour  les  deux  types  de  codes  de  Costas,  avec 
un  maximum  d'environ  -13  dB  k  -14  dB  par  rapport  a  la  valeur 
centrale  de  1' auto-correlation  qui  est  egale  a  la  longueur  du  code. 
Les  resultats  obtenus  pour  la  correlation  croisee  des  codes 
quadratiques  est  un  peu  meilleure  avec  une  moyenne  de  -16  dB. 

Les  codes  etudies  ci-haut  sont  aussi  compares  avec  d' autres 
types  de  codes  qui  ont  deja  fait  l'objet  d' etudes,  par  exemple,  les 
codes  chirp,  step-chirp  et  pseudo-aleatoires  [12,11,1]. 

Nous  avons  compare  le  code  chirp  de  longueur  et  produit  BT 
egaux  k  100  et  les  codes  step-chirp  de  meme  longueur  et  produit  BT 
mais  avec  dix  elements  avec  les  codes  Welch  et  quadratiques  de 
longueur  et  produit  BT  100  ayant  10  ou  100  elements.  Les  codes  de 
Golomb  ne  sont  pas  testes  puisque  qu'ils  n'existent  pas  pour  les 
nombres  d' elements  utilises.  Les  resultats  obtenus  avec  le  code 
chirp  sont  meilleurs  que  les  autres  d'environ  3  dB  a  12  dB.  Les 
resultats  obtenus  avec  le  code  step-chirp  sont  meilleurs  que  ceux 
des  codes  de  Welch  et  quadratiques  d'environ  6  dB  a  15  dB. 

Le  code  chirp  de  longueur  et  produit  BT  dgales  a  121  et  le 
code  step-chirp  de  meme  longueur  et  produit  BT  possedant  11  ele¬ 
ments  ont  ete  compares  aux  codes  de  Golomb  et  quadratiques  de  memes 
longueur  et  produit  BT  ayant  11  Elements.  Les  codes  de  Welch 
n'existent  pas  dans  ce  cas.  Le  code  chirp  est,  ici  aussi,  meilleur 
que  les  deux  autres  avec  un  pic  secondaire  maximum  k  -27.2  dB,  ce 
qui  est  5  dB  de  mieux  que  le  code  de  Golomb  et  15  dB  de  mieux  que 
le  code  quadratique.  Pour  sa  part,  le  code  step-chirp  offrent  des 
resultats  de  8  dB  de  mieux  que  les  codes  de  Golomb  et  de  18  dB  de 
mieux  que  les  codes  quadratiques  avec  un  pic  secondaire  maximum  de 
-30.7  dB. 


Le  tableau  5  montre  les  valeurs  relatives  des  lobes  secon- 
daires  de  la  fonction  d'auto-corr^lation  pour  les  codes  de  Costas 
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(de  Welch  et  de  Golomb) ,  quadratiques,  chirp,  step-chirp  et  pseudo- 
aldatoires.  Ce  tableau  permet  de  comparer  les  resultats  minimums 
et  maximums  obtenus  avec  des  codes  de  longueur  voisine  de  128. 


Code 

Lobes  secondaires 

Minimum 

(longueur; 

61£ments) 

Maximum 

(longueur; 

elements) 

Welch 

-24.0  dB 
(126;  126) 

-13.8  dB 
(126;  18) 

Golomb 

-18.6  dB 
(125;  125) 

-8.9  dB 
(126;  9) 

Quadratique 

-18.0  dB 
(127;  127) 

-0.4  dB 
(128;  128) 

Chirp 

-27.4  dB 
(128;  128) 

-26.3  dB 
(100;  100) 

Step-Chirp 

-30.7  dB 
(121;  11) 

-29.8  dB 
(100;  10) 

Pseudo-aldatoire 

-24.0  dB 
(127;  127) 

-18.0  dB 
(127;  127) 

Tableau  5  Comparaison  de  1' auto-correlation  de  differents  codes. 

Chaque  code  peut  etre  caracterise  par  un  diagramme  d'ambi- 
guite.  Un  diagramme  d'ambiguite  est  une  representation  graphique 
de  1' amplitude  de  la  fonction  d'ambiguite,  laquelle  est  la  reponse 
de  1' auto-correlation  pour  differentes  frequences  Doppler.  Une 
cible  possddant  vine  vitesse  radiale  par  rapport  au  radar  se  trouve 
&  modifier  la  longueur  d'onde  du  signal  reflechi.  La  frequence 
Doppler  est  la  difference  entre  la  frequence  porteuse  du  signal 
transmis  et  celle  du  signal  retourn6,  et  peut  etre  evaluee 
approximativement  par  : 


fd  =  2 •  vd  /  A0  =  2.vd.f0  /  c 

Pour  les  diagrammes  en  forme  de  lame  de  couteau,  comme  ceux 
des  codes  chirp  et  step-chirp  reprdsent6s  h  la  figure  7b,  il  existe 
une  ambiguitd  entre  la  position  r6elle  d'une  cible  en  mouvement  et 
celle  donnde  par  le  corrdlateur.  Pour  les  diagrammes  en  forme  de 
punaise,  comme  ceux  des  codes  binaires,  des  codes  de  Costas  et  des 
codes  h  congruences  quadrat iques  repr6sent6s  ei  la  figure  7a,  il 
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n'est  pas  possible  de  detecter  une  cible  en  mouvement  dont  la  fre¬ 
quence  Doppler  est  superieure  a  1' inverse  de  la  longueur  du  code. 


(a)  (b) 

Figure  7  a)  Fonction  d'ambiguite  d'un  code  de  Costas, 
b)  Fonction  d'ambiguite  d'un  code  chirp. 


Les  diagrammes  d'ambiguite  des  codes  de  modulation  discrete  de 
la  frequence  peuvent  etre  determines  approx imativement  en  super- 
posant  deux  representations  graphiques  de  la  matrice  du  code;  une 
version  demeure  fixe  tandis  que  1' autre  est  deplacee  dans  le  plan 
f requence/temps .  L' amplitude  de  la  fonction  d'ambiguite  est  alors 
proportionnelle  au  nombre  d' elements  qui  coincident  l'un  sur 
1' autre.  Une  cible  en  mouvement  crde  un  decalage  de  la  frequence 
initiale  du  signal  sur  l'axe  des  frequences,  occasionnant  par  le 
fait  meme  un  deplacement  de  la  matrice  sur  l'axe  des  frequences 
Doppler,  tandis  qu'un  decalage  en  distance  implique  un  deplacement 
de  la  matrice  sur  l'axe  temporel. 

La  frequence  des  codes  chirp  et  step-chirp  peut  etre  reprdsen- 
tee  par  une  matrice  avec  une  diagonale.  Un  deplacement  de  la 
matrice  le  long  de  la  diagonale  permet  de  coincider  presque  tous 
les  elements  des  matrices  et  d'obtenir  une  bonne  detection  de  la 
cible.  L' amplitude  de  la  fonction  d'ambiguite  decroit  proportion- 
nellement  avec  ce  deplacement,  lequel  necessite  en  meme  temps  un 
decalage  en  frequence  et  un  decalage  temporel .  Ces  codes  peuvent 
alors  detecter,  avec  une  ambiguite  en  distance,  des  cibles  en 
mouvement  dont  la  frequence  Doppler  est  inf6rieure  A  la  largeur  de 
bande  de  la  modulation  B. 

Les  codes  de  Costas  et  quadratiques  sont  construits  de  fagon 
a  ce  que  la  superposition  d'une  version  decalee  de  la  matrice  sur 
une  version  fixe  ne  permet  de  coincider  que  peu  d' elements.  Ainsi, 
comme  illustre  A  la  figure  7a,  le  diagramme  d'ambiguite  aura  un  pic 
principal  A  l'origine  et  de  faibles  lobes  secondaires  ailleurs  dans 
le  plan  frequence/ temps .  Pour  cette  raison  ces  codes  permettent  de 
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filtrer  les  cibles  en  mouvement  dont  la  frequence  Doppler  est  supe- 
rieure  4  1'  inverse  de  la  longueur  de  1'  impulsion  transxnise. 

•> 

Pour  qu'un  code  puisse  etre  utilise  pour  detecter  une  cible 
dont  la  vitesse  est  super ieure  a  la  limite  imposee  par  le  diagranune  . 

d'ambiguit6,  il  est  ndcessaire  de  modifier  les  coefficients  du  4 

filtre  du  compresseur  de  fagon  h  d6placer  le  diagramme  d'ambigulte 
vers  la  frequence  Doppler  correspondante .  Une  banque  de  plusieurs 
de  ces  filtres  peuvent  Stre  utilises  parallfelement  pour  detecter 
une  plus  grande  gamme  de  vitesses. 


i 


% 


4 


% 
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5.  CONCLUSION 


Ce  rapport  a  permis  d'etudier  deux  differents  types  de  codages 
en  frequence,  les  codes  de  Costas  et  les  codes  a  congruences 
quadratiques .  Ces  codes  ont  ete  evalues  parce  qu'ils  permettent 
d'obtenir  des  diagrammes  d* ambigulte  en  forme  de  punaise  tout  en 
ayant  de  faibles  lobes  secondaires. 

Les  codes  de  Costas  off rent  de  meilleurs  resultats  pour  la 
fonction  d'auto-correiation  du  signal,  surtout  avec  la  methode  de 
construction  de  Welch.  Ceux-ci  sont  en  moyenne  meilleurs  de  5  dB 
par  rapport  aux  codes  de  Golomb  et  de  10  dB  par  rapport  aux  codes 
quadratiques  pour  des  longueurs  voisine  de  128.  Les  codes  quadra¬ 
tiques  ont  par  contre  de  plus  faibles  lobes  pour  la  fonction  de 
correlation  croisee.  Ces  codes  permettent  une  probability  reduite 
de  detection  du  signal  n' ayant  pas  le  merae  code. 

Les  codes  de  Costas  donnent  de  bons  r6sultats  pour  la  hauteur 
des  lobes  secondaires  des  fonctions  d’auto-correiation,  mais  pas 
aussi  bien  que  d'autres  codes  d6je  etudies  comme  les  codes  chirp, 
step-chirp  et  pseudo-aieatoires.  Quant  aux  codes  quadratiques,  ils 
pourraient  etre  appropries  pour  une  utilisation  necessitant  un 
compromis  entre  une  bonne  fonction  d'auto-correiation  et  une  bonne 
correlation  croisee  puisque  les  r6sultats  de  ces  deux  fonctions 
off rent  alors  d'assez  bonnes  possibilites. 
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